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Abstract 

We analyze Optimal Velocity Model (OVM) with explicit delay. The prop- 
erties of congestion and the delay time of car motion are investigated by 
analytical and numerical methods. It is shown that the small explicit delay 
time has almost no effects. In the case of the large explicit delay time, a new 
phase of congestion pattern of OVM seems to appear. 
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I. INTRODUCTION 



In recent years, we proposed a new car-following model "Optimal Velocity Model" (OVM) 
based on a dynamical equation [|l| 

x n (t) = a{V(x n+1 (t) -x n (t)) -x n (t)} , (1.1) 

where t is time and x n is a position of the n-th car. Cars are numbered so that the (n+ l)-th 
car precedes the n-th car. The driver feels the headway x n+ i(t) — x n (t) and determines an 
optimal velocity V(x n+ i(t) —x n (t)). It is best to drive a car with the optimal velocity but in 
general a deviation exists between the optimal velocity and a real one. The driver responds 
to the deviation AV = V{x n+ \(t) —x n (t)) —x n (t) and diminishes it by giving an acceleration 
aAV to the car. The coefficient a expresses the sensitivity of the driver. We call the function 
V "Optimal Velocity Function" (OVF). In previous papers, we have shown how OVM can 
explain behaviors of traffic flow, for example, the transition from a free flow to a congested 
flow, a density-flow relationship, a kind of effective delay of car motion IH-fl- 

On the other hand, the prototype equation of motion of traditional car-following model 

is 

%n — ^o{%n+l ~ %n} > (1-2) 

where Ao is a constant In this model, a driver is thought to react to the stimulus 

proportional to the relative velocity between the previous car and his own car. Equation 
( |1.2|) may be generalized by changing the constant Ao to a function A(x n+ i — x n ) of headway. 
However these models have no physically interesting solution because such equation can be 
integrated easily and be reduced to following equation 

x n = V(x n+ i - x n ) , (1.3) 

where V is a function of headway and V'(x n+ i — x n ) = A(x n+ i — x n ). In car-following 
models, therefore, the introduction of "delay" is necessary and plays an essential role to 
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understand the traffic dynamics @||. Following equation is a typical one which is widely 
used in car-following models. 

x n (t + r) =V(x n+1 (t)-x n (t)), (1.4) 

where r is a delay time of driver's response. The driver senses headway at time t and changes 
the velocity of his car at later time t + r. This delay time r of response has been thought to 
be inevitable because it comes from the driver's physical delay of response to the stimulus 
together with the mechanical response time of a car. In this paper, this r will be called 
"explicit delay time" . 

The notion of explicit delay time r is completely different from that of "delay time of car 
motion" introduced in our previous paper [||. Let us recall the definition of the delay time 
of car motion. Consider a pair of cars, a leader and a follower. Assume the leader changes 
the velocity according to vi = vo(t) and the follower duplicates the leader's velocity but with 
some delay time T, that is, v / = Vo(t — T). Under such a situation we can clearly define the 
delay time of car motion by T. It is known that the observed delay time T of car motion is 
of the order of 1 sec, but the known physical or mechanical response time r is of the order 
of 0.1 sec. In the previous paper we confirmed that the equation ( |1.1| ) really produces T of 
order 1 sec. 

We clarified that OVM can describe the properties of traffic flows or the behaviors of 
cars fairly well without any explicit delay time r. However there exists the delay time of 
response of driver for a fact. The explicit delay time r should be included in the dynamical 
equation in order to construct realistic models of traffic flow. It is a natural question what 
kind of effect appears in the traffic flow or in the car motion if we introduce the explicit 
delay in the equation(|i~T|). 

In this paper we investigate the following equation 

x n (t + r) = a{V(x n+1 (t) -£„(*)) -x n (t)} . (1.5) 

In order to our analysis be more concrete, we use the parameter a = 2.0 (1/sec) and the 
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function V which are phenomenologically determined in previous papers by the observed 
data on Japanese motorways [|T^-[12|. 

V(Ax) = 16.8 [tanh 0.0860 (Ax -25) + 0.913] , (1.6) 

in which the unit of length and time are 'meter' and 'second' respectively. 

The plan of this paper is as follows. In section 2 we discuss the global properties of traffic 
flow in OVM with the explicit delay. In section 3 we investigate more microscopic property, 
that is, the delay time of car motion. First we discuss within a linear approximation and 
next evaluate the delay times of car motion in various cases by numerical simulations. As a 
special case, the car motion under the traffic signal is also treated. In section 4 we show the 
new feature of OVM with the explicit delay. The final section is devoted to summary and 
discussion. 

II. PROPERTY OF TRAFFIC FLOW IN OVM WITH EXPLICIT DELAY 

A. Linear Analysis 

In this section we investigate OVM with the explicit delay time r of driver's response 
described by equation ( |1.5|) . 

First we analyze the linear stability of iV-car system on a circular lane of length L. 
Obviously, the homogeneous flow solution of equation (|1.5| ) is given by 

x£\t) = V(b)t + nb, b = L/N . (2.1) 

To see whether the solution (|2.1| ) is stable or not, we add a small perturbation 

x n (t)=x^(t)+y n (t) . (2.2) 

From equation Q2.2p and equation (|1.5[), we can calculate a linearized equation with respect 
to y n (t) 

y n {t + T) = a{fAy n {t)-y n {t)}, (2.3) 
3 



where / = V'(b) and Ay n = y n +i — y n - A complete set of solutions is given by 

Vjn(t) = exp(i<x,n + iuijt), (2.4) 
where otj = 2-nj/N for j — 1, 2, 3, • • - , N and Uj satisfies the equation 
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In equation fl2.5|) , variables are combined to be dimensionless. The condition that each 
solution yj n {t) becomes marginally stable is Im Uj = 0. For convenience of explanation, we 
will omit the mode-index j and treat a as a continuous variable. The condition Im u = 
gives 'critical curves' for each ar in (f/a, a) plane, where f/a is a radial coordinate and a 
is an angular coordinate. Mode solutions yj n (t) are represented by a point (f/a,aj) on a 
circle //a=const. 

Three critical curves for ar = 0, 0.2 and 0.4 are shown in figure 0, in which a reference 
circle represents mode solutions for f/a = 0.75. The modes staying outside (right-hand 
side) of the critical curve are unstable. Figure [I] shows that a homogeneous flow state with 
a parameter f/a = 0.75 is an unstable state. From figure [I], it is found that unstable modes 
increase as the explicit delay time r becomes large. This situation looks similar to a case 
the sensitivity a becomes small in the original OVM [|IJ. Though there seems to be some 



relationship between the sensitivity a and the explicit delay time r as indicated in Ref. |L3 
an analytical relation has not been clarified yet. 



B. Numerical Simulations 

The effect of the explicit delay in the congestion formation can be evaluated by numerical 
simulations. In previous papers ]1],0, we investigated the property of traffic flows in a circuit. 
It is found that when the car density exceeds a critical value, a homogeneous traffic flow 
becomes unstable and makes a phase transition to a congested flow. After enough time, 
the congested flow becomes stationary and shows an alternating pattern of high density 
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(congestion cluster) and low density regions. Each velocity and headway inside high (low) 
density regions always take common values which are determined only by the sensitivity a 
and OVF independently of any other conditions. The motion of each car can be shown in a 
'phase space' (Ax, x), and the trajectories draw a single "hysteresis loop", a kind of limit 
cycle. Figure [2] shows typical hysteresis loops for sensitivity a = 2.0 and 2.8. Two turning 
points C = (Ax c ,v c ) and F = (Axf,Vf) correspond to the high and low density region for 
a = 2.0 and C and F 1 for a = 2.8. We found the congestion pattern moves backward on 
the circuit with a constant velocity (vfAx c — v c Ax f) / (Ax f — Ax c ), which is given by the 
intersection of x axis and the line connecting two turning points C and F. Therefore the 
property of such congested flows is almost decided by two points C and F of hysteresis loop. 

From numerical simulations, we recognize no qualitative difference in the behavior of the 
traffic flow between the cases with and without the explicit delay, if r is not so large. Figure 
|3] shows hysteresis loops for r = 0, 0.1 and 0.2, that is, ar = 0, 0.2 and 0.4. The changes of 
hysteresis loops are similar to those for the case that the sensitivity a becomes small in the 
original OVM []]]. Therefore it seems that the explicit delay time r, which is not so large, 
does not play any essential role in the congestion formation. In other words, the effect of 
the explicit delay can be almost compensated by the change of sensitivity a. 

Obviously, this is not the case for a very large ar. Figure |4] shows examples for ar = 
0.6, 0.8 where critical curves are inside the referenced circle f/a— 0.75. In original OVM 
instability always comes from long range modes (a ~ 0), that is, short range modes (a ~ 7r) 
are always stable. In the case ar > 0.6, however, there exist various cases in which all 
modes become unstable or short range modes only become unstable. In such cases, the 
instability starts from all modes or from short range modes. It is interesting to see what 
kind of phenomena emerge in such cases. An example shall be discussed in section 4. 



III. DELAY TIME OF CAR MOTION 
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A. Linear Analysis 



In this section, we investigate the delay of car motion in order to see the effect of the 
explicit delay from a more microscopic viewpoint. First, we analyze the delay of car motion 
with the linear approximation. 

Consider a pair of a leader and it's follower where the leader moves with the velocity 
v(t) and the follower replicates the motion of the leader after the time interval T, that is, 
the follower's velocity is given by v(t — T). In this case we can define the delay time of car 
motion as T. 

Let the position of the leader at time t be y(t) and that of it's follower x(t) which obeys 
equation (|1.5|) , that is, 

x(t + r) = a{V{y(t) - x(t)) - x(t)} . (3.1) 

Starting from the situation with headway b and velocity V(b), 

y (t) = V(b)t + b, x (t) = V(b)t, (3.2) 

we investigate the response of the follower to a small change X(t) of the leader: 

y(t)=y {t) + X{t), x(t) = x (t)+£{t). (3.3) 

Inserting above equations into equation ( |3 . 1|) and taking a linear approximation, we get 

i(t + r) + a£(t) + af^it) = af\(t), (3.4) 

where / = V'(b) is again a derivative of OVF at headway b. If one takes X(t) = e tuJt , the 
solution is given by 

= 1 ■ • If 1 iun 2/ f ^ ■ ( 3 - 5 ) 

1 + tu/j — e %u)T uj z /aj 

This is rewritten as 

e(t) = \i\ e^- T) , (3.6) 
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where 



_ 1 _-, auj — oj 1 sin(a;r) 
T = — tan — - — 

10 a j — ui l cos(uit) 



lei 



,CL> r> ,ijP . , U) . ,UP »n 

1 + - 2 - 2 — cos cur + -smcur + — 2 
f af f af 



(3.7) 
(3.8) 



First let us consider the case |u;| is sufficiently small (no/ f <S 1, o>/a <g 1). It will be discussed 
later whether this condition is satisfied or not in the realistic situation used in equation ( |1.6|) . 
Then we have 



iei = i, t 



r 



(3.9) 



Here we take the general expression of X(t) which is expressed as follows. 



X(t)= / X(uj)e luJt du; . 



(3.10) 

X(u>) is assumed to be nonzero only for uj small enough. Then we find the follower's response 
becomes 



f (t) = / X(u)e iuj ^ T) du = X(t - T) 



(3.11) 



that is, 



x(t) = V(b) + £(t) = V(b) + X(t -T)= y{t - T), 



(3.12) 



with T of equation (p. 9|). 

As a result we conclude that for sufficiently slow and small change of leader's velocity, 
the delay time T of motion of the follower becomes 1/ / (the inverse of derivative of OVF at 
corresponding headway), independently of the explicit delay time r of driver's response. 



B. Simulations for Homogeneous Flows 



Next we will carry out numerical simulations to investigate the effect of the explicit delay 
in homogeneous traffic flows. The validity of the conditions u a, f can be checked also. 

7 



We make simulations starting from homogeneous flows with various headways and add a 
small disturbance to the first car. The delay time of car motion is estimated between 10th 
car and 11th car when the disturbance propagates there. 

In table |, we summarize the results of numerical simulation. In the cases where the 
homogeneous flow is stable, the delay time T of car motion is almost equal to 1/f and the 
explicit delay has no effect. The cases Ax = 20, 25, 30 correspond to the unstable situation. 
The measurement of the delay time T is carried out before the disturbance becomes large. 
The results show that the assumption uj <C a, / is not valid. Even in such cases the explicit 
delay does not affect T. 

C. Simulation for Congested Flows 

In this subsection, we treat the car motion in a stationary congested flow, where linear 
analysis is no more valid obviously. In the previous paper [|J we have shown that the delay 
time T of car motion is the inverse of the gradient of line which connect two turning points 
(C and F in figure |2|). This is a natural extension of the statement obtained by the linear 
analysis: "The delay time of car motion is the inverse of derivative of OVF at corresponding 
headway" . 

Our task here is to carry out similar numerical simulations with the explicit delay. After 
the congestion pattern becomes stable, all cars behave in the same manner expressed in 
figure [| We can estimate the delay time T from the time interval of the motion of successive 
two cars, which is equivalent to the gradient of line connecting two turning points of the 
hysteresis loop. Table [Ij shows the results of simulations for r = 0, 0.1, 0.2. 

The table clearly shows that the change of T is rather small compared to the change of 
t. Therefore the main contribution of the delay of car motion comes from the structure of 
OVM itself and not from the explicit delay. The r-dependence of T appears only through 
the change of turning points of the hysteresis loop. In other words, the effect of the explicit 
delay is similar to the change of the sensitivity a and is not essential in the same as the 
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previous section. 



D. Simulations for Car Motion under a Traffic Signal 

In this subsection we study the delay of motion of cars starting from a traffic signal. 
Though this may be a special case compared to previous subsections, experiments to observe 
the delay time have often been done in this situation. 

Numerical simulations are carried out as follows. First a traffic signal is red and all cars 
are waiting with headway 7 (m), at which OVF (|1.6|) becomes zero. At time t — 0, the 
signal changes to green and cars start. 

Figures |5] and ^| show the velocities of several cars in a queue for the cases of r = and 
t = 0.2 (sec), respectively. It can be seen that cars with large car number (7th or more) 
behave almost in the same manner as its preceding car. 

We can estimate the delay time T from the behavior of velocities of 7-10th cars. Table 
[TTI| shows the delay time of car motion for various r. Again we find that the delay time T 
has a small dependence on the explicit delay time r. To see whether this is general or not, 
we carried out another simulations with the initial headway 3 (m). For this purpose, OVF 
( |1.6|) is changed to take zero for Ax < 7 (m). We show the results in the third column of 



table |n| which again show obviously a small dependence of T on r. 

Hitherto we concerned the definition of delay time of car motion given in the section 1: 
if velocities of two successive cars are given by v(t) and v(t — T) respectively, the delay time 
of car motion is T. This definition is valid only for the case that the motions of two cars are 
similar. As is seen from figures |5] and [| the first several cars move in the different manner, 
because the headway of the first car is infinite but that of other cars are relatively small. 
In order to explore the delay time of car motion in such case we will propose an another 
definition. For example, we can define the delay time as the interval between the time when 
the preceding car starts and the time when the next car starts. Though there are many 
other possibilities, the above definition looks rather natural. 
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Figure [7| shows the delay time of car motion by the new definition. Obviously the data 
approach to a certain values as the car number becomes large. The limits of the delay 
times in this definition are the same values as those in the previous definition. It should be 
mentioned that the explicit delay time r is simply added to the delay time T of car motion 
for first a few cars. This effect dissipates after several cars start. 

From these results, we can conclude that the explicit delay time r contributes directly 
to the delay time of car motion only for such a restricted for the motion of first a 

few cars starting from the traffic signal. In general case, the contribution of r to T is rather 
small and is similar to the contribution from the change of the sensitivity a. 

IV. NEW FEATURES OF OVM WITH EXPLICIT DELAY 

In this section we show new features which exist only in OVM with the explicit delay. 

A. Overshoot Phenomenon 

We investigated the motion of cars controlled by a traffic signal in the previous section. 
For small r the motions of cars are not so different from those for r = 0. For large r, 
however, we can see a transitional overshoot of velocity, that is, a excess and a gradual 
decrease of velocity. As a typical case, the motions of cars for r = 0.3 are shown in figure || 
We have carried out many numerical simulations by changing r and found that the overshoot 
phenomenon begins at r = 0.19 (sec). 

B. Upper Bound of r 

First we note that the explicit delay time r is understood as the summarized effect 
coming from delays of physical and mechanical response. Therefore too large value will not 
be permitted from observations. There exists, however, more restrictive bound, which has 
a origin in the equation of motion (|1 . 1|) of OVM. 
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We consider a homogeneous equation of the linearized equation (|3.4|) in the leader- 
follower system: 

Z{t + T)+a4{t)+af£{t) = 0. (4.1) 

£(t) gives a perturbative motion of the follower when the leader moves in a constant velocity. 
In order that the two body system is stable, £ (t) must vanish as time develops. We see that 
£(t) = e lujt is a solution of equation ( |4.1|) , with u satisfying 

- cu 2 e iujT + iau + af = 0. (4.2) 

The marginally stable condition Imw = becomes 

en = Ksin(«;) , fr = Kcot(ft) , (4.3) 

where k = Re u; r. By eliminating k, we can find the upper bound of r for given a and /. 
Though we could not solve equation (|4.3j ) analytically, the upper bound r m is found to be a 
monotonic decreasing function of both a and /. 

The value of be evaluated numerically. For the sensitivity a = 2.0 (1/sec) and the 



maximum value of / = 1.44 (1/sec), which is read off from OVF ( |1.6|) , the corresponding 
upper bound r m is 0.44 (sec). If r > r m (sec) in OVM with above a and /, the car cannot 
follow the constant velocity motion of the leader. Thus r m should be understood as the 
upper bound of the explicit delay time in order that OVM is meaningful as a model of 
traffic flow. 



C. New Congestion Pattern 

Inside the above upper bound of the explicit delay time, some curious phenomena emerge 
in traffic flow as the explicit delay time becomes large. If such phenomena should be regarded 
as unrealistic, the upper bound will be taken at a smaller value. 

Figure (|9|) shows a snap shot of headway after enough simulation time. The conditions 
of the simulation are as follows: total car number N = 100, circuit length L = 2500 (m) 
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and explicit delay time r = 0.22 (sec). There we can see small congestion clusters or rapid 
change of velocity between 15th and 60th car. This pattern looks like a intermediate pattern 
before the congestion is formed completely. However, in contrast to the case of r = 0.20 
(sec) where such a pattern or small congestions disappear as time goes, the pattern has very 
long life and may never disappear in the case of r = 0.22 (sec). This pattern occupies a 
larger region as r increase. 

Next we take r to be a larger value 0.4 (sec). Figure [K| shows the hysteresis loops for 
t = and 0.4. Here we note that OVF ( |1.6|) takes negative value continuously for Ax < 7 
(m) and therefore cars can move backward (without collisions). Because such behaviors 
of vehicles are obviously unrealistic, it seems natural to set the upper bound of r to the 
transition point at which this hysteresis loop appears. 

As shown in figure pi], the profiles of hysteresis loops are qualitatively different. Moreover 
the hysteresis loop for a = 2.8 is larger than that for a = 2.0 in the case of r = 0.4 in contrast 
to the case of r = 0. We also note that the relaxation time for r = 0.4 is of the order of 
10 3 ~ 10 4 times that for r < 0.2. The differences of hysteresis loops and relaxation times 
seem to suggest an existence of a new phase. However, there exists another possibility: the 
stationary state indicated by this hysteresis loop is artificial due to finite size effects and a 
new phase does not exist. The congestion pattern changes continuously around r ~ 0.22 
(sec) and we cannot find a definite transition point. It is a future work whether this pattern 
indicates the existence of new phase or not. 

V. SUMMARY AND DISCUSSION 

In this paper we investigated the properties of Optimal Velocity Model with the explicit 
delay of driver's response. The effects of the explicit delay are very small, if the delay time 
is small: r < 0.2 (sec). The effects are similar to the change (reduction) of the sensitivity a, 
and therefore the explicit delay does not play an essential role. This fact should be compared 
to the traditional car-following models, in which the delay of driver's response has played a 
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significant role. The equation of motion of traditional car-following model becomes trivial, 
if the delay time is zero. 

For large explicit delay time r, the traffic flow behaves in a different manner. If r < 0.2 
(sec), the properties of congestion clusters are similar to that for r = 0. For r > 0.2 (sec), 
the stationary pattern of the traffic flow does not consist of only such congestion clusters 
but confused patterns. For r > 0.3 (sec), the traffic flow becomes stationary but congestion 
clusters are never formed. 

In OVM, there is an upper bound of the explicit delay time, which comes from the 
condition that the equation of motion is meaningful. The upper bound, however, becomes 
small, if we require the existence of stable congestion clusters. 

From this work, we can obtain an indication on a phenomenological study. In this paper 
we clarified the notions of the delay time r of driver's response and the delay time T of 
car motion. However the meaning of the delay time of response and its effect are model- 
dependent. In traditional car-following models, the delay time r seems to be merely a fitting 
parameter and so we can take any value for r. Moreover, the delay time often takes different 
value in each term. In OVM, the delay time r is not free and the observed value decided 
by experiments will give a criterion whether OVM with OVF ( |1.6| ) is valid or not. Here we 
note that the contribution of the delay of driver's response to the delay of car motion is 
very small. The delay of car motion, therefore, has it's root just in the dynamical equation 
itself. This fact suggests the difficulty to determine the delay time r of driver's response by 
measuring the delay time T of car motion. Therefore r must be measured directly by other 
experiments. 
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FIG. 1. Solid line, dashed line and dotted line show critical curves for ar = 0, 0.2 and 0.4, 

respectively. A circle of diamond marks represents mode solutions for f /a = 0.75. 
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FIG. 2. "Hysteresis loops" for a = 2.0 and a = 2.8. Each line connects two turning points of 
each hysteresis loop. A tanh-type curve represents OVF (|1.G|). 
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FIG. 3. Hysteresis loops for r = 0, 0.1 and 0.2. A tanh-type curve represents OVF (|1.6[). 
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FIG. 4. Solid line and dashed line show critical curves for ar = 0.6 and ar = 0.8, respectively. 

A circle of diamond marks represents mode solutions for f /a = 0.75. 
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FIG. 5. Motions of No. 1-11 cars for r = 0. Each curve shows the velocity of each car. 
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FIG. 6. Motions of No. 1-11 cars for r = 0.2. Each curve shows the velocity of each car. 
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FIG. 7. Solid line connects delay times (time intervals) of 2-1 1th cars for r = 0. Dashed, dotted 



and dashed dotted lines connect those for r = 0.1, 0.2 and 0.3, respectively. 



40 I 1 1 1 1 1 r 




2 4 6 8 10 12 14 

time (s) 



FIG. 8. Motions of No. 1-11 cars for r = 0.3. Each curve shows the velocity of each car. 
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FIG. 9. A snapshot of velocities at t = 10000 (sec). Diamond marks represent velocities of cars. 
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FIG. 10. Hysteresis loops for r = and r = 0.4. As a reference, two cases, a = 2.0 and a = 2.8, 
are shown. A tanh-type curve represents OVF ( |1.6| ). 
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TABLES 



A y> ( y-v-i j 

L — \Ju I 111 \ 


f _1 fsl 
J W 


- t T=0 1 & J 


- t T=0.1 V & / 


- L T=0.2 


10 

-LU 


9 6497 




9 fi 


9 fi 


1 ^ 






1 .OO 


1 3^ 
l.OO 


90 


&9£9 


u.yo 


u.yo 


n 


25 


0.6921 


0.85 


0.87 


0.89 


30 


0.8282 


0.95 


0.95 


0.95 


35 


1.3434 


1.35 


1.35 


1.35 


40 


2.6427 


2.6 


2.6 


2.6 


50 


13.101 


13 


13 


13 



TABLE I. Delay times of car motions in homogeneous flows. 



T simulation 



0.0 0.94 

0.1 0.96 

0.2 0.99 

TABLE II. Delay times of car motions in congested flows. 



r (s) T for headway 7m (s) T for headway 3m (s) 
0.0 1.10 1.26 
0.1 1.10 1.26 
0.2 1.11 1.25 
0.3 1.12 1.26 



TABLE III. Delay times of car motions in queues starting from a traffic signal. 
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